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Abstract. The chemistry of the early Universe is a fascinating field of study. Even in the 

absence of any elements heavier than lithium, a surprising degree of chemical complexity proves 

" ' to be possible, giving the topic considerable interest in its own right. In addition, the fact that 

I , molecular hydrogen plays a key role in the formation of the first stars and galaxies means 

f^ ■ that if we want to understand the formation of these objects, we must first develop a good 

rvj ' understanding of the chemical evolution of the gas. In this review, I first give a brief introduction 

. : to the chemistry occurring in the gas prior to the formation of the first stars and galaxies, 

^ and then go on to discuss in more detail the main chemical processes occurring during the 

^~~i gravitational collapse of gas from intergalactic to protostellar densities, and how these processes 

f — . influence the final outcome of the collapse. 
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T 1. Introduction 

O I Chemistry in the early Universe has proved to be an interesting topic of study for 

J,^ ■ several reasons. Before the first stars and galaxies formed, the elemental composition 

\^ 1 of the gas was very simple - it consisted only of hydrogen and helium (and their iso- 

[/3 ' topes) plus a tiny trace of lithium - and yet this is sufficient to give rise to a sur- 

C^ prising degree of chemical complexity. In addition, the interaction between the gas 

and the cosmic microwave background (CMB) plays an important role in regulating 

y—i . the chemistry of the gas, and moreover one which may lead to observable changes in 

> ; the CMB (see e g. iMaoh. Melchiorri fc Tostilll994HDubrovich. Baikova fc Khaikinil2008 : 

''^ . ISchleicher et al.ll2008 ). Furthermore, understanding the chemical evolution of the gas at 

^_ ' these early epochs is vital if one wants to understand the formation of the first stars and 

fy\ . galaxies, given the crucial role played by molecular hydrogen (H2) in this process. 

There is a wide range of topics that one could discuss in a review of the chemistry of 
the early universe, ranging from the physics and chemistry of the recombination epoch 
at redshift z > 800 to the reionization of the Universe at 2: < 15. However, to keep this 
review to a manageable length, I will focus on only a few key issues, and will not discuss 
either the recombination epoch or the process of cosmologi c al reio nization. Both of these 
topics are treated in detail elsewhere (see e.g. iFendt et al.l ( 2009f) and references therein 



o 



s 



p% ' for a lengthy discussion of th e physics of the recombination epoch, or the recent review 



by [Morales fc Wvithd ( 20101 ) for a discussion of the epoch of reionization) 



In the next Section, I review what wc know about the chemical evolution of the IGM 
from the end of the recombination epoch until the time at which the earliest protogalaxies 
form. I follow this up in Section |3] by discussing the chemical evolution of the gas within 
these early protogalaxies. I show how an understanding of the gas chemistry can allow us 
to determine the mass scale of the first star-forming minihalos, outline the main chemical 
changes that occur in the gas during its collapse to protostellar densities, and highlight 
the important role played by molecular hydrogen throughout this process. 
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2. Pregalactic chemistry 

2.1. The chemical composition oj the gas at the end of the recombination epoch 

At the end of the recombination epoch, at a redshift z ^ 800, the chemical makeup of the 
gas in the Universe - what we might refer to as pregalactic gas - was very simple. The 
main constituent of the gas was neutral atomic hydrogen, accounting for roughly 75% of 
the total mass. A small amount of hydrogen remained ionized, wit h the ratio by number 
of ion ized to neutral hydrogen being roughly H"*"/!! ~ 2.5 x 10~^ (jWong. Moss fc Scott 
2008t note that this value continues to decline at lower redshift, decreasing by about an 
order of magnitude by the time we reach z = 100). Neutral atomic helium (primarily ^He, 
but also including a small fraction of '^He) accounted for most of the remaining mass, but 
in this case the abundance of singly or doubly ionized helium was very close to zero, as 
helium recomb ination comes to an end at a significantly higher redshift than hydrogen 
recombination ( Wong. Moss fc Scottll2008l ) . Deuteri um was also present, with a fra ctional 
abundance relative to hydrogen D/H ~ 2.5 x 10~^ (jCvburt. Fields fc Olivell2008l ). and a 
ratio of D"*" to D that was essentially the same as the ratio of ioni zed to neutral hydrogen. 



Final ly, there was a trace amount of lithium (Li/H ~ 5 x 10~^°; ICvburt. Fields fc Olive 
2008[) . the majority of which was in the form of Li+. The fractional abu ndances of all other 
ionized or molecular species were smaller than 10^^^ at this point (JAlizadeh fc Hirata 



2.2. The chemical evolution of the IGM 

Starting from this simple beginning, many other chemical species were produced as the 
Universe evolved. Detailed chemical models for the evolution of pregalactic gas have 
been constructed by a number of authors (see e.g. Stancil. Leop & Dalgarnqll996l 1998 : 
Galli fc Pallall998l : (Lepp. Standi, fc Dalgarnol2002HPuv fc Sign orc 2002; Schleic her et al, 
20081 ). but in this discussion I will focus on a few general principles that can help us to 
understand the chemical evolution of the gas at these early times. 

As a starting point, it is useful to draw a distinction between what I will call "primary" 
species, i.e. species that are formed directly from chemical reactions amongst the main 
chemical constituents of the gas (H, H"*", D, D"*", He and Li"*"), and "secondary" species 
that are formed mainly from reactions involving one or more of the primary species. 
Good examples of primary species arc H^ , formed via the reaction 



H + e ^H +7, 



(2.1) 



or H2 , formed via the reaction 



H + H+^H++7. (2.2) 

The prime example of a secondary species is H2 , as direct formation of H2 via the reaction 

H + H^H2+7 (2.3) 



occurs at a negligible rate (JGould fc Salpeteiill963l ). and most H2 in metal- free gas forms 
via the reactions 



H- + H ^ H2 + c- 



(2.4) 



and 



H++H^H2+H+. (2.5) 

The evolution of the fractional abundance of a primary species X is largely governed 
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by the relative sizes of three important timescales: the formation time, defined as 



Rf, 



'orin,X 



(2.6) 



where nx is the number density of species X and -Rform.x is the formation rate per unit 
volume of X; the photodissociation time, 

ipd.x = -f^, (2.7) 

-Kpd,X 

where i?pd,x is the photodissociation rate of species X; and the Hubble time 

tu = J^y (2.8) 

where H{z) is the Hubble constant. 

At high redshift, photodissociation of most molecules and molecular ions by the CMB 
occurs rapidly, and the photodissociation time is small. Similarly, the high gas density 
means that the formation timescale of most species is also small. In this regime, both 
timescales are typically much shorter than the Hubble time, and hence the abundance 
of species X evolves until it reaches an equilibrium value set by the balance between 
formation and photodissociation (in which case iform.x — ^pd.x)- Since photodissociation 
is very effective at high redshift, this equilibrium value is generally very small. 

As we move to lower redshifts, the cosmological background density decreases, since 
p{z) oc (1 + z)'^, and hence the formation time of any given primary species tends to 
increase. However, the energy density of the CMB decreases even more rapidly, prad oc 
(1 + z)^, and the CMB temperature decreases as Tcmb oc (1 + z). What this means 
is that we eventually reach a redshift at which the number density of photons capable 
of photodissociating species X starts to fall off exponentially, as the threshold energy 
for photodissociation starts to fall within the Wien tail of the CMB spectrum[|J Once we 
reach this point, ipd rapidly increases, and we are left with a situation in which iform is by 
far the shortest of the three timescales. The abundance of species X therefore increases 
rapidly until tform,x — ^h- Beyond this point, the abundance of X does not increase 
significantly, as the time required would be longer than the age of the Universe, and we 
speak of the abundance of the species "freezing out" at some asymptotic value. 

The evolution with redshift of secondary species such as H2 is harder to generalize in 
this way, as it often depends on the behaviour of more than one primary species, each of 
which passes through the three stages noted above at a different redshift. For example, 
H2 formation via the H^ ion becomes possible at a redshift z ~ 400, but H2 formation 
via the H~ ion is suppressed until z ^ 100, owing to the fact that H~ has a much lower 
binding energy than H^ (0.755 eV for the former, 2.65 eV for the latter), and hence 
can be destroyed by photons with significantly smaller energies. The H2 abundance in 
the pregalactic gas therefore passes through two periods of relatively rapid growth, one 
occurring around z ^ 400 and the second around z ^ 100. Following this second period 
of growth, the fractional abundance of H2 is of the order of 10~^, and H2 is the most 
abundant molecule in the gas. However, this is still far smaller than the typical amount 
of H2 produced within the earliest star-forming protogalaxies, and hence although the 
chemical evolution of the pregalactic gas is interesting, it is of limited influence: almost 

f This simple picture is complicated somewhat by the presence of non-thermal photons pro- 
duced during hydrogen recombination that slightly d istort the spectral shape of the CMB 
(jSwitzer fc Hiratall2005l : iHirata fc PadmanabharJl2006D . but even in this case the general prin- 
ciple is the same. 
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all of the chemistry affecting the formation of the first stars and galaxies occurs within 
the protogalaxies themselves, as I explore in more detail in the next section. 



3. Protogalactic chemistry 

As we move to lower redshifts, the distribution of matter in the Universe becomes 
increasingly inhomogeneous. The small perturbations in the dark matter distribution 
that were present during the recombination epoch grow with time, evolving initially as 
S oc t^'"^, but eventually entering the regime in which their evolution becomes non- linear. 
These perturbations undergo runaway gravitational collapse and virialization, forming 
dense, gravitationally bound structures that are often referred to as dark matter halos. 
In the ACDM model, the earliest halos to form do so on the smallest scales, with larger, 
more massive halos forming at later times. 

Very little gas accumulates within the earliest halos, as it is prevented from undergoing 
gravitational collapse by its own thermal pressure, and also by the significant relativ e 
velocity that exists between the gas and the dark matter rTseliakh ovich fc Hiratall201Cll ). 
However, once the dark matter halos become sufficiently massive, their gravitational 
attraction dominates over these effects and substantial quantities of gas flow into them, 
being heated by adiabatic compression and shocks in the process. In regions where the 
relative velocity between gas and dark matter is very s mall, this occurs once t he mass 



of the dark matter halo exceeds M^i^ ~ 2 x 10'' Mq (|Naoz fc Barkanal 120071 ). but in 



more typical regions where there is a relative velocity between gas and dark matter 
of several kilometers per second, the m inimum halo mass is higher, Mmin ~ 10^ M, 



( Tseliakhovich. Barkana fc Hiratal 120101 ). Low- mass halos with masses close to Af, 







often referred to as "minihalos" to distinguish them from their considerably more massive, 
low-redshift brethren - are therefore the first objects that have the potential to form stars. 
Whether or not the gas within them can actually form stars depends on how rapidly it 
can cool. This in turn depends on the chemical evolution of the gas within the minihalos, 
a subject that I address below. 

3.1. H2 formation and the onset of cooling 

As gas falls into a typical minihalo, it is heated up to a temperature close to the virial 
temperature of the minihalo 

(where this expression assumes tha t the gas is predominantly atomic hydrogen, and that 
z ^ 10; see lBarkana fc Loebll20Ql l for a more general expression). In the post-shock gas. 



the electron fraction decreases due to radiative recombination, but at the same time II2 
forms, primarily via reactions 12 . II and 12 .41 The amount of II2 that forms i n the gas can be 
estimated using a simple toy model for the chemistry first introduced bv lTegmark et al 



(1997 ). We start by assuming that radiative recombination is the only process affecting 

the electron abundance, and writing the rate of change of the electron number density 

as 

drio ,„ -, 

—— = -/erocno«H+, (3.2) 

at 

where n^ is the number density of electrons, n^+ is the number density of protons, and fciec 

is the recombination coefficient. If we assume that ionized hydrogen is the only source 

of free electrons, implying that Uc = nH+ , and that the temperature remains roughly 
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constant during the evolution of the gas, then we can solve for the time evolution of the 
electron fraction: 

X = ; , {'J-J) 

1 + KrccntXo 

where x = ric/n, n is the number density of hydrogen nuclei, and xq is the initial value 
of x. If we now assume that all of the H2 forms via the pair of reactions given above, and 
that the only process competing with reaction 12.41 for the H~ is mutual neutralization 

H++H"^H + H, (3.4) 

then we can write the time evolution of the H2 fraction, xhj = ^Ha /"-j ^s 

—rf- = ^raa^nnPad, (3.5) 

where A:ra is the rate coefficient of reaction l2.1[ and Pad is the probability that any given 
H~ ion will be destroyed by associative detachment (reaction 12. 4|) rather than by mutual 
neutralization. Given our assumptions above, this can be written as 

Pad = T —, > (3-6) 



where k^d is the rate coefficient for reaction 12.41 and fcmn is the rate coefficient for reac- 
tion 13.41 If we again assume that Ue = nH+ , and in addition assume that the H2 fraction 
and the fractional ionization are small, so that riu — n, then the expression for pad can 
be simplified to 



"^ad 

Substituting this into Eauation l3.5[ we obtain 



Pad = ( 1 + ^X ) . (3.7) 



kj-i^xnn ( 1 + -^^x ] . (3.8) 



dt \ fcad 

We can solve for the time evolution of the H2 fraction 

Kra , / 1 + ^Jofcnm/fSad + ^ I i-ccc \ /„ n\ 

^H2 = 1 — In T— ] 77 , (3.9) 

"Tcc \ -L + 2:0 "rnn / "^ad / 

where tree = (^rec'^2;o)~^ is the recombination time, and where we have once again 
assumed that tih ~ n- If the initial fraction ionization is small, so that Xq <C ka^d/kmm 
then we can write this expression in the simpler form 

XH, = ^ In (1 + f Aroc) . (3.10) 

This analysis tells us two important things. First, it shows us that the H2 fraction 
evolves only logarithmically with time: most of the H2 forms within the first few recom- 
bination times, with H2 formation at later times being strongly suppressed owing to the 
loss of the free electrons from the gas due to recombination. Second, it shows us that 
the amount of H2 that forms depends on the ratio of the rate coefficients for H" forma- 
tion (reaction 12. ip and H"*" recombination, fcra/fcroc- The value of this ratio is strongly 
temperature dependent, and can be written approximately as 

^ c. 10-«Ti -^ (3.11) 

The amount of H2 that forms in the gas therefore increases significantly with increasing 
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tempera ture. As the H2 coo ling rate is also a strongly increasing function of temperature 
(see e.g. iFlower et al.lbood ). one finds that there is a reasonably sharp divide occurring 



at a critical virial temperature Tcrit '~ 1000 K between minihalos that contain gas that 
can cool rapidly and form stars, and minihalos that contain g as that barely cools over a 
Hubble time (see e.g. lTegmark et al.lll997tlYoshida et al.ll2003 ). If we convert this critical 



virial temperature into a corresponding critical minihalo mass using Equation l3.11 we find 
that 

Merit ^ 7.6 X lO^/i-i (^) M0. (3.12) 

At z > 40, this value is smaller than the minimum mass Mmin ~ 10^ Mq required 
to overcome the effects of the streaming of the gas relative to the dark matter, which 
therefore sets the minimum mass scale for a star-forming minihalo. At z < 40, however, 
Merit > Mmin, and it is the amount of H2 that can form within the gas that is responsible 
for setting the minimum mass scale. One implication of this result is that there will be 
a population of small minihalos that forming at z < 40 that contain significant gas 
fractions, but that do not form stars, because their gas is unable to cool in less than a 
Hubble time. These small starless mini halos may be im portant sinks for ionizing photons 
during the epoch of reionization ([Haiman. Abel fc Madau 2001i) . 

3.2. The initial collapse phase 

In minihalos where H2 cooling is efficient, the gas cools and collapses under its own self- 
gravity, collecting at the centre of the minihalo. The evolution of the gas at this stage 
depends on the amount of H2 it is able to form, which in turn depends upon the initial 
ionization state of the gas. 

During the formation of the very first Population HI stars, the initial fractional ion- 
ization of the gas is the same as the residual ionization in the intergalactic medium, i.e. 
xo ~ 2 X 10"** at z ~ 30. In this case, the amount of H2 that forms in the gas is typically 
enough to cool it to a temperature of T ^ 200 K but not below. At this temperature, 
deuterated hydrogen, HD, forms easily from H2 via 

H2+D+^HD + H+, (3.13) 

but its destruction via the inverse reaction 

HD + H+^H2+D+ (3.14) 

is starting to become ineffective, owing to the fact that this reaction is endothermic by 
around 0.04 eV. At T = 200 K, chemical fractionation therefore boosts the HD/H2 ratio 
by a factor of 20 compared to the cosmic deuterium-to-hydrogen ratio. This boost to 
the HD abundance, together with the fact that HD is a m uch more effective coolant 
than H2 at low temperatures (see e.g. iFlower et al.l I2OOOI) , allows HD to become an 



important coolant at 200 K, but detailed studies have shown that the extra cooling 
provided by HD is not enough to reduce the temperature significantly below 200 K 



(jBromm. Coppi fc Larsonll2002f ). Therefore, H2 continues to dominate the cooling and 
control the further evolution of the gas. In this scenario, the collapse of the gas is greatly 
slowed once its density reaches a value of around 10** cm~^, corresponding to the critical 
density ncrit, at which the rotational and vibrational level populations of H2 approach 
their local thermodynamic equilibrium (LTE) values. At densities higher than this criti- 
cal density, the H2 cooling rate per unit volume scales only linearly with n (compared to 
a quadratic dependence, Ahj oc n^ at lower densities), while processes such as compres- 
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sional heating continue to increase more rapidly with n. As a result, the gas temperature 
begins to increase once the density exceeds ricrit- 



Gas reaching this point in the collapse enters what iBromm. Coppi fc LarsonI (|2002i ) 
term a "loitering" phase, during which cold gas accumulates in the centre of the halo 
but only slowly increases its density. This loitering phase ends once the m ass of cold gas 
that has ac cumulated exceeds the local value of the Bonnor-Ebe rt mass (|Bonnoiill956 : 
Ebertlll955l) . given in this case by ( Abel. Brvan. fc Normanll2002 ) 

(3.15) 



MBE=i40T3/^n-i/2MQ, 



which for n - 10'* cm'^ and T ~ 200 K yields Mbe ~ 1000 M©. Once the mass of cold 
gas exceeds Mbe, its collapse speeds up again, and becomes largely decoupled from the 
larger-scale behaviour of the gas. The next notable event to occur in the gas is the onset 
of three-body H2 formation, which is discussed in Section 13.31 below. 

If the initial fractional ionization of the gas is significantly higher than the residual 
fraction in the IGM, then a slightly different chain of events occurs. A larger initial 
fractional ionization implies that tree is shorter and hence leads to a logarithmic in- 
crease in the amount of H2 formed after a given physical time. The higher H2 fraction 
produced in this way allows the gas to cool to a slightly lower temperature, thereby 
boosting the HD abundance to a significantly higher level than is reached in the low 
ionization case. This allows HD to take over as the dominant coolant, driving the temper- 
atu re down further. A number of s tudies have examined this HD-dominated reginie (see 
e.g. Nakamura fc Umemura 2002 : Naeakura fc Omukail 2005t Johnson &: Brommll2006 



Yoshida et al.ll2007l: iMcGreer fc Brvanll2008t iKreckel et al.ll2010[) , and have shown that 
the minimum temperature reached in this case can be as low as the CMB temperature, 
and that the higher critical density of HD, ricrit.HD ~ 10^ cm~'^, means that the gas does 
not reach the loitering phase until much later in its collapse. Once the gas does reach this 
phase, however, its subsequent evolution is very similar to that in the low-ionization case 
discussed above. Cold gas accumulates at n ~ ncrit until its mass exceeds the Bonnor- 
Ebert mass, which in this case is Mbe ^ 40 M© if T = 100 K and n = 10^ cni^^. Once 
the gas mass exceeds Mbe, the collapse speeds up again, and the gas begins to heat up. 
Aside from the substantial difference in the size of A^be, the main difference between the 
evolution of the gas in this case and in the low ionization case lies in the fact that in 
the high ionization case, the gas reheats from T r^ 200 K or below to T ~ 1000 K much 
more rapidly than in the low ionization case. As we shall see later, this period of rapid 
heating may exert an important influence on the ability of the gas to fragment. 

Several different scenarios have been identified that lead to an enhanced fractional 
ionization in the gas, and that potentially allow the gas to reach the HD-dominated 
regime. Gas within minihalos with Tvir > 9000 K will become hot enough for coUi- 
sional ionization of hydrogen to supply the necessary electrons, but halos of this size 
will typical l y hav e at least one star- forming progenitor (1 Johnson. Greif fc Bromml 12008 : 



GreifetaL 



20101 ). and hence typically will have already been enriched with metals. 
Other possibilities include the formation of protogalaxies in "fossil" HII regions, i.e. 
regions th at are recombining after having been ionized by a previous Popul ation HI 
protostar ( Oh fc Haimanll2003t iNagakura fc Omukail 12005 : Yoshida et al.l 2007 ). or that 
are irra diated with a significant flux of X-rays (e.g. Glo ver fc Brand! 12003 ) or cosmic 
ravs (iStacv fc Bromm 2007 : Jasche. Ciardi fc Ensslin 20071 ). However, recent work by 
Wolcott-Green fc HaimanI (|201l[ ) has shown that the HD forming within a primordial 
minihalo is very susceptible to p hotodissociation by even a weak extragalactic UV back- 
ground (see also lYoshida. Omukai fc Hornguist 2007 ). and so it is possible that gas will 
reach the HD-dominated regime in only a few minihalos. 
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3.3. Three-body H2 formation 

Once the collapsing gas reaches a density of around 10^-10^ cm^'', its chemical makeup 
starts to change significantly. The r eason for this is that at these den sities, the formation 
of H2 via the three-body reactions (jPalla. Salpeter. fc Stahleiill983[ ) 



H + H + H 


^H2 


+ H, 


I + H + H2 


^H2 


+ H2 


I + H + He 


^H2 


+ He 



(3.16) 
(3.17) 
(3.18) 

starts to become significant. These reactions quickly convert most of the hydrogen in 
the gas into H2. At the same time, however, they generate a substantial amount of 
thermal energy: every time an H2 molecule forms via one of these three-body reactions, 
its binding energy of 4.48 eV is converted into heat, and so at these densities, three-body 
H2 formation heating can become the main process responsible for heating the gas. For 
this reason, even though the abundance of H2 , the dominant coolant during this phase 
of the collapse, increases by more than two orders of magnitude, the gas typically does 
not cool significantly. Indeed, the temperature often increases. 

One major uncertainty in current treatments of the gas chemistry at this point in the 
collapse is exactly how quickly the gas becomes molecular. Although reaction l3.16l is the 
dominant source of H2 at these densities, the rate coefficient for this reaction is poorly 
known, with published values differing by almost two orders of magnitude at 10 00 K, 
and by an even larger factor at lower temperatures (Glover 2008 : Turk et al.ll201ll) . The 
effects of this uncertainty have recently been studied bv I Turk et al.l (2011). They show 
that it has little effect on the density profile of the gas, and only a limited effect on the 
temperature profile. However, it has far more significant effects on the morphology of 
the gas and on its velocity structure. Simulations in which a high value was used for the 
three-body rate coefficient find that gravitational collapse occurs more rapidly, and show 
that the molecular gas develops a much more flattened, filamentary structure. Significant 
differences are also apparent in the infall velocities and the degree of rotational support. 
The effect that this has on the ability of the gas to fragment (see Section 13.51 below) is 
not yet understood. 

3.4. The final stages of collapse 

As the gas collapses to even higher densities, several important events occur which affect 
its ability to cool. The first of these occurs at a density of around 10^° cm~^, where 
the main rotational and vibrational lines of H2 start to become optically thick. This 
reduces the effectiveness of H2 cooling, lea ding to an ongoing ris e in the gas temper 



ature. In one-dimensional simulations (e.g. lOmukai fc Nishil 119981 : lOmukai et al.l 11998 ; 



Ripamonti et al.l I2OO2I : iRipamonti fc Abell 12004 ). it is possible to treat optically thick 
H2 cooling accurately by solving the full radiative transfer problem. These models show 
that although the optical depth of the gas becomes large at frequencies corresponding to 
the centers of the main H2 emission lines, the low continuum opacity of the gas allows 
photons to continue to escape through the wings of the lines, with the result that the 
H2 cooling rate is suppressed far less rapidly as the collapse proceeds than one might 
at first expect fsee lOmukai et al.lll998l for a detailed discussion of this point). In three- 
dim ensional simulations , more approximate treatments of the cooling are necessary (see 
e.g. lYoshida et al] I2OO6I: iTurk. Alpel fc O'SheallJOOOt IClark et al.l 1201 lal) . but the gen- 
eral findings are the same - H2 cooling is suppressed to some extent, but still provides 
enough dissipation of energy to allow the collapse to continue without a large increase in 
the temperature. 
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A second important event occurs once the number density reaches n ^ 10 cm~ . 
At this density, colhsion-induced e mission from H2 ( the inverse of the more famihar 
colhsion-induced absorption; see e.g. lFrommholdlll993 ) becomes important. An isolated 
H2 molecule has no dipole moment, but the activated complex formed during the collision 
of two H2 molecules can have a dipole moment, and hence can emit radiation via dipole 
transitions. In principle, this can occur in gas of any density, but the probability of a 
photon being emitted in any given collision is very small, owing to the short lifetime of 
the collision state (A^ < 10~ ^^s at the temperatures relevant for Pop. Ill star formation; 



Ripamonti fc Abelll2004[) . For this reason, collision-induced emission (CIE) becomes 



see 

an important process only at very high gas densities. Another consequence of the short 
lifetime of the collision state is that the individual lines associated with the dipole transi- 
tions become so broadened that they actually merge into a continuum. This is important, 
as it means that the high opacity of the gas in the rovibrational lines of II2 does not 
significantly reduce the amount of energy that can be radiated away by CIE. There- 
fore, once the gas reaches a sufficiently hig h density, CIE becomes the dominant form of 



cooling, as pointed out b y several authors (jOmukai fc Nishilll998t iRipamonti et al.ll2Q02 
Ripamonti fc Abelll2004l ). 



The most detailed s tudy of the effects of CIE cooling on the collapse of primordial 
gas was carried out by IRipamonti fc Abell (J200J). They showed that CIE cooling could 
actually become strong enough to trigger a thermal instability. However, the growth rate 
of this instability is longer than the gravitational free-fall time, meaning that it is unlikely 
that this process can drive fragmentation during the initial collapse of the gas. 

The phase of the collapse dominated by CIE cooling lasts for only a relatively short 
period of time. The g a.s becomes optically thick in the continuum o nce it reaches a 
density n ~ 10^^ cm"^ (jOmukai fc Nishi|[l998l: IRipamonti fc Abelll2004[) . which strongly 
suppresses any further radiative cooling. Once this occurs, the gas temperature rises until 
it reaches a point at which the H2 begins to dissociate. At these densities, this occurs at 
a temperature T ~ 3000 K. The dissociation of H2 slows the temperature rise for a while, 
as most of the energy released during the collapse goes into dissociating the H2 rather 
than raising the temperature. However, once the H2 is exhausted, the temperature once 
again begins to climb steeply. The thermal pressure in the interior of the collapsing core 
rises rapidly and eventually becomes strong enough to halt the collapse. At the point at 
which this occurs, the size of the dense core is around 0.1 AU, its mass is around 0. 01 MfT, 



and its mean density is of order 10^" cm ^ ( Yoshida. Omukai fc HernquistI l2008[ ) . It is 



bounded by a strong accretion shock. This pressure-supported, shock-bounded core is 
the structure that we identify as a Population HI protostar. 

3.5. Formation and fragmentation of the protostellar accretion disk 

Most numerical simulations of the formation of Population HI stars stop at or before 
the point at which the protostar forms, owing to the increasingly small timesteps that 
need to be taken to follow the evolution of the gas. A s significant fragmentation o f the 
gas during the initial collapse is uncommon (see e.g. iTurk. Abel fc O'Sheal l2009l who 
find that fragmentation into a binary system occurs in only one out of five cases) , most 
studies have assumed that the gas surrounding the protostar does not fragment at later 
times either, but instead is simply accreted by the existing protostar, either directly or 
via a pr otostellar accretion d isk. A good overview of the results of these studies can be 
found in lBromm fc Yoshidal (201lJ). 

However, in the past few years, a number of simulations have been performed that 
attempt to directly model the evolution of the gas after the formation of the first Popula- 
tion HI protostar using a technique developed for studies of contemporary star formation 
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(seee .g.lClarket alJl2008t IStacv. Greif fc Broinml[2010t IClark et alJ[2Qllal lbl: iGreif et al. 



20111 : ISmith et al.ll201iri . Gravitationally bound regions of gas that become smaller than 
some pre-selected size scale are replaced by what are usually termed sink particles (in 
an SPH code; see e.g. iBate. Bonnell fc Pricelll995l) or sink cells (in a grid code; see e.g. 



Federrath et al.l 120101 ). These sinks can accrete gas from their surroundings and con- 



tinue to interact gravitationally with the surrounding gas, but allow one to neglect the 
very small-scale hydrodynamical flows that would otherwise force one to take very small 
numerical timesteps owing to the Courant condition. 

One of the main results of these simulations is the demonstration that the first protostar 
to form does not simply accrete all of the available gas. The infalling gas forms an 
accretion disk around the protostar, but this disk is unable to transfer gas inwards, onto 
the protostar, rapidly enough to balance the rate at which new material is falling onto 
the disk. The disk therefore grows in mass until it becomes gr avit at ionally unst ab le , at 
which point it begins to fragment, forming multiple protostars ( Clark et al.ll2011bl ). The 
gravitational interaction between these protostars and the surrounding gas can trigger 
further fragmentation, and close encounters between protostars may lead to mergers, or 



in so me cases to protostars being ejected from the center of the minihalo (jGreif et al 
20111) . 



Molecular hydrogen plays a crucial role in this process. The midplane temperature of 
the protostellar accretion disk is typically 1500-2000 K, and the disk is maint ained at this 



temperature primarily by a combination of H2 line cooling and CIE cooling (|Clark et al 
l2011b) . If the molecular hydrogen were not present, the main disk coolant would be 
bound- free emission from H~ , which becomes effective only at a temperature T ^ 6000 K. 
The disk would therefore be considerably hotter, allowing it to transfer gas more rapidly 
onto the central protostar, and making it less susceptible to gravitational instability. 
In this case, previous stu dies have shown that the disk would probably remain stable 
and w ould not fragment (JTan fc McKeel 120041: iTan fc Blackmanlbooi iMaver fc Duschll 
200Rah . 



As well as regulating the disk temperature, molecular hydrogen also provides the 
cooling necessary to allow the fragments forming in the disk to survive as distinct, 
gravitationally-bound objects. Studies of the fragmentation of self-gravitating accretion 
disks have shown that in order fo r fragments to survive, they must be able to cool within 
roughly half of an orbital period (JGammiell2001l ). Fragments that cannot cool this quickly 
undergo a thermal bounce and are then sheared apart by disk rotation, rather than con- 
tinuing to collapse. In the case of the fragments forming in a Pop. Ill accretion disk, 
two key processes provide the necessary coolin g: CIE cooling at n < 10^^ cm~^ and the 



collisional dissociation of II2 at n > 10^® cm ^ (jClark et al.ll2011bl ) 



Finally, there are some indications that the chemical state of the gas at early times 
during its collapse exerts a substantial influence on the degree of fragmentation that 
occurs at much later times. As we have already seen, if HD cooling becomes efficient 
during the collapse of the gas, then the gas can cool to a significantly lower temperature 
than in the standard H2-dominated case. However, a consequence of this is that the gas 
undergoes a more rapid period of reheating once the HD molecules reach LTE, during 
which it has a very stiff effective equation of state. This leads to a pronounced loss of 
small-scale structure from the gas, which appears to reduc e the amoun t of fragment ation 
that occurs at later times in the evolution of the minihalo ( Clark et aLluOllai: iGreif et al.l 
201 ll) . 
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